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This article presents experimental results on properties of waves propagating in an unbounded 
and a bounded photonic crystal consisting of metallic cylinders which are arranged in a triangular 
lattice. First, we present transmission measurements of plane waves traversing a photonic crystal. 
The experiments are performed in the vicinity of a Dirac point, i.e., an isolated conical singularity of 
the photonic band structure. There, the transmission shows a pseudodiffusive 1/L dependence, with 
L being the thickness of the crystal, a phenomenon also observed in graphene. Second, eigenmode 
intensity distributions measured in a microwave analog of a relativistic Dirac billiard, a rectangular 
microwave billiard that contains a photonic crystal, are discussed. Close to the Dirac point states 
have been detected which are localized at the straight edge of the photonic crystal corresponding to 
a zigzag edge in graphene. 

PACS numbers: 42.70.Qs, 73.22.Pr, 42. 25. Fx 



I. INTRODUCTION 

The experimental realization of graphene 1 , a one- 
atom-thick sheet of carbon atoms arranged in a honey- 
comb crystal lattice 2,3 , attracted a lot of interest due to 
its peculiar electronic band structure which exhibits a lin- 
ear dispersion at low energies. There the band structure 
closely resembles the spectrum of a relativistic massless 
spin one-half particle 4 . As a consequence, even though 
the Fermi velocity of electrons in graphene is only c/300 
with c the speed of light, the electronic excitations are 
governed by the Dirac equation. The peculiar electronic 
properties lead to interesting physical phenomenal. 

The linear dispersion relation is entirely due to the 
threefold rotational symmetry of the hexagonal graphene 
lattice 6 . Therefore, graphene with its Dirac spectrum is 
not an exception. In general, systems comprising waves 
that propagate in a spatially periodic potential with the 
same symmetry properties can posses a Dirac spectrum. 
One example are so-called photonic crystals 7 , optical 
analogs of the common ones. It was shown theoretically 
in Ref£ that under certain conditions photonic crystals 
with a triangular lattice geometry, depicted schematically 
in Fig. [lj can exhibit the linear Dirac dispersion rela- 
tion. The dependence of the frequencies / of propagating 
electromagnetic waves on the components (k x ,k y ) of the 
quasi-momentum vector exhibits a band structure which 
is similar to that for the electronic energies in graphene. 
The left hand side of Fig. [2] shows the first two bands 
obtained for a two-dimensional photonic crystal consist- 
ing of metallic cylinders of radius R = 0.25a, where a is 
the lattice constant. At the cylinders, the electric field 
strength vanishes, that is, obeys the Dirichlet boundary 
condition. The band structure was calculated by numer- 
ically solving the scalar Helmholtz equation for the elec- 
tric field strength using the finite difference method de- 
scribed in Ref. 9 . The two bands touch each other at the 
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FIG. 1: (Color online) Schematic view of the triangular lattice 
of metallic cylinders (black) . The voids between the cylinders 
are marked by colored circles. For graphene the red (bright) 
and the blue (dark) voids correspond to the two atoms which 
generate the two independent triangular sublattices. The ar- 
rows indicate the two directions TM and TK considered in the 
transmission experiments. The radius of the cylinders equals 
R = 0.25 a, where a is the lattice constant. The dashed lines 
indicate the positions of the metallic walls used in the ex- 
periments described in Sec. Ill Bl to realize a microwave Dirac 
billiard. 



corners of the hexagonal Brillouin zone depicted in Fig. [2] 
where they have the shape of cones. Thus, there the fre- 
quencies of propagating modes depend linearly on the 
distance of the quasi momentum from the touch points. 
The Brillouin zone comprises two non-equivalent corners 
K and K', so-called K points, that cannot be mapped 
onto each other by basis vectors of the reciprocal lattice. 
The band diagram closely resembles that of the electronic 
energies in graphene 10 . Accordingly, as in grapheneii the 
touch points are referred to as Dirac points. 

The photonic crystal also exhibits a hexagonal config- 
uration. In fact, three neighboring cylinders forming a 
triangular cell of the photonic crystal lattice shown in 
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Fig. [I] host quasibound states^ 2 - localized at their center. 
The cell thus can be considered as an open resonator and 
the triangular lattice can be regarded as being built of 
coupled resonators, the voids between the cylinders act- 
ing as the atoms forming the two independent triangular 
sublattices that generate the honeycomb structure. In 
Fig. [I] the voids corresponding to the two sublattices are 
marked with blue and red circles, respectively. The void 
structure terminates with a zigzag edge along the TK and 
with an armchair edge along the TM direction, respec- 
tively. These are the two most common types of bound- 
aries arising from the hexagonal structure of graphene. 

In the vicinity of the Dirac points the Schrodinger 
equation for the electrons in graphene ^ 13 i 14 and the 
Helmholtz equation governing the electromagnetic wave 
propagation in the photonic crystal 8 reduce to two copies 
of the Dirac equation for massless relativistic spin one- 
half particles, one copy holding around K and the other 
around K'. Thus, photonic crystals can be used to in- 
vestigate properties of graphene, as was done in ana- 
log experiments with microwave photonic crystals^ - — 
and also with sonic crystals^. There phenomena di- 
rectly related to the linear dispersion relation around 
a Dirac point were investigated. We consider here two 
of them, the extremal transmission and the formation of 
edge state s 13 ! 19 ! 21 . Both are properties of the waves prop- 
agating through a photonic crystal and crucially depend 
on the structure of the crystal edges. 

The Dirac dispersion relation affects the transport 
properties of photonic crystals in a remarkable way. It 
was shown theoretically in Refs j 19 i 2Q that near the Dirac 
point the intensity of electromagnetic waves transmitted 
through a photonic crystal is inversely proportional to 
the thickness of the crystal. This is different from the 



ballistic behavior and the exponential decay observed, 
respectively, in a band region and a gap region, but 
rather reminiscent of the diffusion behavior of waves trav- 
eling through a disordered medium. A similar behavior 
was predicted for the conductivity of an ideal graphene 
strip 21,22 . We present below results on the experimen- 
tal investigation of the transmission of electromagnetic 
waves traversing a microwave photonic crystal. The mea- 
surements (see Sec. Ill A|) were performed close to a Dirac 
point, which was located as described in Ref. 17 . 

In further experiments (see Sec. Ill Bp we investigated 
edge effects ^ 23 1 24 on the properties of the wave func- 
tions 25,26 of a microwave analog of Dirac billiards. These 
experiments were motivated by the results obtained for 
the spectral properties of graphene flakes, also called 
graphene billiards 24 , which have been shown to depend 
on the geometry of their edges 25 27 . For a graphene flake 
terminated with armchair edges the density of states van- 
ishes at the Dirac point such that a gap opens up in the 
band structure, whereas zigzag edges yield a nonvanish- 
ing density. Indeed, in the vicinity of the Dirac point 
graphene flakes which terminate with zigzag edges or 
with both armchair and zigzag edges exhibit edge states, 
which are localized close to the zigzag parts of its bound- 
ary- 2 ^""—. This property is attributed to the boundary 
condition s 24 ! 31 which differ for both types of edges. To 
check the existence of such states in photonic crystals of 
finite size, one with the structure shown in Fig. [1] was en- 
closed by a metallic frame and wave functions were mea- 
sured in the vicinity of a Dirac point. While in graphene 
flakes electrons are confined by a potential barrier, the 
walls of the metallic box prevent microwaves from leav- 
ing the photonic crystal and correspond to an infinitely 
high potential barrier in the quantum analog. 
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FIG. 2: (Color online) Left: band structure f(k x ,k y ) of the 
triangular lattice shown in Fig. [1] as function of the quasi- 
momentum components (k x , k y ). It was computed by numer- 
ically solving the scalar Helmholtz equation for the electric 
field strength imposing the Dirichlet boundary condition at 
the metallic cylinders. Shown are the first and the second 
band. Blue (darkest) color corresponds to the lowest, red 
(brightest) color to the largest values of / within a band. 
The hexagonal Brillouin zone with the two non- equivalent K 
points K and K' is indicated. Right: band structure com- 
puted numerically by varying the quasi-momentum vector 
along the path TMKr defining the boundary of the irreducible 
Brillouin zone depicted together with the Brillouin zone in the 
inset. Adopted from Ref.—. 



II. EXPERIMENTS WITH THE MICROWAVE 
PHOTONIC CRYSTAL 

Analog experiments with flat, cylindrical microwave 
cavities have been performed to model generic proper- 
ties of quantum chaotic systems and of chaotic quan- 
tum scattering systems 3 - 2 -"— . This is possible in a range 
/ < /max of the excitation frequency /, where the wave 
length is longer than twice the height of the cavity. There 
the electric field vector E is perpendicular to the bottom 
and the top plates of the cavity and obeys the Dirichlet 
boundary condition at the side walls of the cavity. Ac- 
cordingly, the vectorial Helmholtz equation reduces to a 
scalar one which is mathematically equivalent to the two- 
dimensional stationary Schrodinger equation for a non- 
relativistic quantum billiard of corresponding shape— - —. 
Due to this analogy flat microwave cavities are commonly 
called microwave billiards. We used microwave photonic 
crystals consisting of metallic cylinders which were ar- 
ranged in a triangular lattice (see Fig. [1]) and squeezed 
between two metal plates to model single particle prop- 
erties of graphene. At the cylinders the electric field 
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strength vanishes, i.e., it obeys the Dirichlet boundary 
condition. Thus, for frequencies / < / max the setup cor- 
responds to a quantum multiple-scattering system. The 
experiments are described in Sec. Ill Al and Sec. IIIB| re- 
spectively. 



A. Extremal transmission at the Dirac frequency 

For the experimental investigation of the transmission 
properties of a photonic crystal with a Dirac spectrum 
we used the setup shown in Fig. [U The photonic crystal 
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FIG. 3: Experimental setup (not to scale). The photonic crys- 
tal consisted of 960 copper cylinders arranged in a triangular 
lattice. A microwave horn antenna with an opening angle 
of 2A# = 28° was used as emitting antenna and the signal 
transmitted through the crystal was received by a waveguide 
antenna. The crystal was placed symmetrically with respect 
to the antennas and oriented such that normal incidence of the 
emitted waves coincided with the TK direction, in distinction 
to the experiments described in Ref.— where it corresponded 
to the TM direction. 

consisted of 960 cylinders with radius R = 5 mm and 
height h = 8 mm. To ensure a good electrical contact 
and thus reproducibility of the measurement each cylin- 
der was screwed to the top and the bottom plate. The 
resulting photonic crystal contained 40 rows of 24 cylin- 
ders and had side lengths of, respectively, 400 mm and 
859 mm. The height of the cylinders defined the maxi- 
mal frequency /max = 18.75 GHz up to which the system 
was two dimensional. A two-dimensional horn antenna 
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FIG. 4: Examples of transmission spectra measured with the 
setup shown in Fig. [3] in TK direction through a lattice with, 
respectively, 8, 18, and 36 rows of cylinders. The dashed line 
marks the Dirac frequency. 



was used as emitting antenna and a waveguide antenna 
received waves transmitted through the crystal. Both 
were tightly screwed to the top plate and attached to 
waveguide-to-coaxial adapters. The transmission spec- 
tra were obtained with an Agilent PNA-L N5230A vec- 
torial network analyzer (VNA). It measures the ratio of 
the power P5 of the signal received at the waveguide an- 
tenna and the power P a of that emitted from the horn 
antenna and their relative phase. The modulus square 
of the scattering matrix element Sb a from antenna a to 
antenna b is given as 



\s ba (f)\ 2 = ^- 



(1) 



The horn antenna emits a beam with divergence 2A#, 
where A6 = 14° with respect to the direction indicated 
by the arrow in Fig. [3j which impinges the crystal situ- 
ated in its near field. The orientation of the crystal was 
chosen such that normal incidence coincided with the TK 
direction (see Fig. [2]). In contrast, in the experiments de- 
scribed in Ref. 17 it corresponded to the TM direction. 

The transmission spectra shown in Fig. [4] were taken 
for different configurations of the experimental setup. 
First, the transmission spectrum was measured for the 
setup with 40 rows and then, step by step one row was 
removed from each armchair edge of the crystal. The 
three panels show exemplarily the spectra measured for 
crystals consisting of 8, 18, and 36 rows of cylinders, re- 
spectively. In all spectra the transmission is suppressed 
below 11.8 GHz. There, the band structure exhibits a 
gap (see right panel of Fig.[2j), and electromagnetic waves 
cannot propagate in the photonic crystal. Furthermore, 
the transmission spectra show a characteristic dip around 
the Dirac frequency fp = 13.797 GHz, which was deter- 
mined experimentally in Ref. 17 and is marked in Fig. |4] 
with a dashed line. 

To reveal the 1/L scaling, with L denoting the crystal 
thickness, of the transmitted power \Sb a (f)\ 2 predicted 
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FIG. 5: The product of the transmission \Sb a (f)\ 2 through 
the photonic crystal and its thickness L in units of the lattice 
constant a as function of L for different frequencies, from top 
to bottom for 13.18 GHz, 13.797 GHz, and 18.7 GHz. The 
filled circles are experimental data. They are joined with 
dashed lines to guide the eye. The solid lines result from 
the fits of a linear, a constant and an exponential function, 
respectively, to the data points. 



near the Dirac point, we show in Fig. [5] the quantity 
(L/a)\Sba(f)\ 2 versus L/a for three different frequencies. 
In all three panels the filled circles are the experimental 
data. They are joined by dashed lines to guide the eye. 
The solid curves depict the average behavior. The top 
panel shows (L/a)\Sb a \ 2 versus L/a for / = 13.18 GHz, 
where the transmission is nonvanishing. In this frequency 
regime L\Sb a \ 2 grows linearly with L, that is, the trans- 
mission |5& a | 2 itself does not depend on L, as expected in 
band regions. The oscillations around the mean are at- 
tributed to standing waves arising between the antennas 
and the interfaces of the photonic crystal. The bottom 
panel of Fig. [5] shows the variation of (L/a)\Sb a \ 2 with 
L/a for / = 18.7 GHz. This frequency lies within a stop 
band. As expected the transmission decreases exponen- 
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FIG. 6: Double- logarithmic plot of the transmission |Sb a (/)| 2 
through the photonic crystal as function of L at the Dirac 
frequency fn = 13.797 GHz. The filled circles are experimen- 
tal data. The solid line shows the prediction for the average 
behavior obtained by evaluating Eqs. (4)- (6) of Ref. 20 . The 
fit of Eq. J2} describing the oscillations 18 yields the dashed 
line. 



tially with the crystal thickness L. The middle panel 
shows (L/a)\Sba\ 2 versus L/a at the Dirac frequency 
f D = 13.797 GHz. The average behavior of L\S ba \ 2 
is constant over a wide range of the crystal thickness. 
Thus, the transmission of plane waves excited with the 
Dirac frequency and traveling in directions close to the 
TK direction, i.e., with \k\\\ « 0, where k\\ is the wave 
vector component parallel to the interface of the crystal, 
scales with \/L. This is similar to the diffusion behav- 
ior of waves propagating through a random medium. As 
was pointed out in Ref^ 9 - this phenomenon of extremal 
transmission is due to the linear Dirac dispersion rela- 
tion. Indeed, in the vicinity of a Dirac point the local 
density of states vanishes only exactly at the Dirac fre- 
quency, whereas in a band gap it is zero for a finite range 
of excitation frequencies. As a consequence, the observ- 
ability of the 1/L scaling crucially depends on the orien- 
tation of the photonic crystal with respect to the incident 
plane waves. When emitting plane waves excited with the 
Dirac frequency in TM direction, as was done in previous 
experiments^^, the band structure has a partial band 
gap (see right panel of Fig. [2|). Accordingly, the trans- 
mission \Sba\ 2 decreases exponentially with the crystal 
thickness L, that is, the decay length of the evanescent 
waves is finite. In contrast, it is infinite when exciting the 
Dirac point, e.g., by sending plane waves at the Dirac fre- 
quency in TK direction onto the crystal. 

The 1/L scaling near a Dirac point has been predicted 
for photonic crystals in Ref. 19 . For its revealment the 
solutions of the Helmholtz equation in free space were 
matched at the interface of the crystal to those of the 
Dirac equation, which is applicable inside the photonic 
crystal in the vicinity of the Dirac point. Furthermore, 
conservation of |fc|||, that is, Snell's law, and of the flux 
through the boundary had to be taken into account. In 
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Fig. [6] we show the data of the middle panel of Fig. [5] in a 
double-logarithmic plot together with the prediction for 
their average behavior (full line) obtained by evaluating 
Eqs. (4)- (6) of Ref. 20 . Here, the range of angles ±A# 
of the propagation direction with respect to the TK di- 
rection of the waves impinging the crystal is AO = 14°. 
This yields for the range ±A of the displacement of the 
wave vector k from the K point at K, that is, of \k\\\, the 
value A = A67rf D /c ~ (7v/3)a~ 1 . For log 10 (L/a) > 0.7. 
Hence, since the thickness of a layer equals a/2 (see 
Fig. [1]), at least 10 layers are needed, in order that the 
average transmission | S& a | 2 through the crystal shows the 
1/L scaling predicted for L/a > l/(aA) = 3/tt. It has 
been observed before in two-dimensional sonic crystals^ 
and in photonic crystals consisting of dielectric rods 15 . 
There, in addition, oscillations of the transmission were 
found for an arrangement of the crystal with respect to 
the wave source as in Fig. [3] and attributed to the scatter- 
ing between the valleys at K and K'. Oscillations are also 
visible in the double-logarithmic plot of the transmission 
shown in Fig. [6] They are, however, not as pronounced 
as those presented in Refs j 15 i 18 . In distinction to these 
experiments, successively two rows, instead of one row of 
cylinders were removed for the measurement of the data 
presented in Figs. [5] and [6] Nevertheless, the period of the 
oscillations is well reproduced by the equation 15 (dashed 
line) 

\S ba \ 2 = Aa/L (1 - exp[-LA]) ([1+g + g cos[|K|L - *]) , 

(2) 

where A, g and $ are fitting parameters. The fit of 
Eq. (j2j) to the date yielded A = 0.015, g = -0.55 and 
$ « tt/3. 



B. Edge states in a rectangular Dirac billiard 

A motivation for the experiments presented in this sec- 
tion were the results obtained for the spectral proper- 
ties of graphene flakes^22~— and Ref4£. There Berry 
and Mondragon considered the Dirac equation for mass- 
less spin one-half particles confined to a two-dimensional 
domain by introducing an infinite-mass term along its 
boundary, i.e., a so-called neutrino billiard. Close to 
the Dirac frequency the propagation of electromagnetic 
waves in a photonic crystal can be effectively described 
by the Dirac equation. Thus, enclosing a microwave pho- 
tonic crystal by a metallic frame, where the electric field 
strength obeys the Dirichlet boundary condition and the 
outgoing current vanishes, corresponds to introducing an 
infinite-mass term in the associated Dirac equation. Ac- 
cordingly, below / = /max this provides an experimen- 
tal setup for the investigation of properties of relativistic 
quantum particles put inside an infinitely high potential 
barrier. The reflection of electromagnetic waves at the 
walls of the metallic box leads to an excitation of states 
in both the valleys at K and K' of the associated band 
structure. Their coupling depends on the edge struc- 



FIG. 7: (Color online) Photograph of the Dirac billiard with 
the top plate removed. The photonic crystal with a triangular 
lattice geomatrey is enclosed by a rectangular brass frame 
with side lengths of, respectively, 218 mm and 420 mm. It 
consisted of 273 copper cylinders of radius R = 5 mm and 
height h = 8 mm, had a lattice constant a = 4R and was 
squeezed between two copper plates. The two white crosses 
mark the positions of the antennas. 



ture. Therefore, in distinction to neutrino billiards, such 
systems are governed by two Dirac equations, which are 
coupled depending on the shape of the metallic box en- 
closing the crysta l 13 ! 26 . They are generally called Dirac 
billiard o 25 i 27 . Figure [71 shows a photograph of the exper- 
imental realization of a rectangular Dirac billiard, where 
the top plate has been removed. The photonic crystal 
had the same lattice constant a = AR and consisted of 
273 copper cylinders with the same radius R = 5 mm 
and height h = 8 mm as in the scattering experiment de- 
scribed in Sec. Ill Al It is surrounded by a brass frame with 
side lengths of, respectively, 218 mm and 420 mm. In the 
experiments all cylinders and the frame were squeezed 
with screws between the top and bottom plates to ensure 
an optimal electrical contact. To measure transmission 
spectra wire antennas of 1 mm diameter were lowered 
into the resonator through 3 mm wide drillings in the 
top plate. They reached only 2 mm into the interior of 
the cavity to minimize the disturbance of the electric field 
modes excited there. The positions of the antennas are 
marked by white crosses in Fig. 

The top panel of Fig. [8] shows a measured transmis- 
sion spectrum of the Dirac billiard. The positions of 
the resonance maxima yield the eigenvalues of the Dirac 
billiard. The stop bands below 11.8 GHz and between 
18 and 20 GHz are well pronounced. As in the scatter- 
ing experiment described in Sec. Ill Al we observe a dip 
around the Dirac frequency fo = 13.797 GHz, which 
is marked with an arrow. The bottom panel of Fig. [8] 
shows a magnified section of the transmission spectrum 
around the Dirac point. The spectrum contains sharp 
resonances and many overlapping ones. The latter ren- 
der the extraction of the complete sequence of resonance 
frequencies impossible. However, complete sequences of 
the resonance frequencies are needed for the study of 
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FIG. 8: Top panel: transmission spectrum measured between 
two wire antennas. Bottom panel: magnification of the spec- 
trum in a region around the Dirac frequency. The arrows 
mark the resonance frequencies for which the measured elec- 
tric field intensity distributions are shown on Fig. [9] The 
labels of the resonances refer to the panels in Fig. [9] 



their fluctuation properties. This corroborates the ne- 
cessity of the measurement of the spectra at supercon- 
ducting condition ^ 41 ! 42 in order to reduce dissipation of 
microwave power in the walls of the metallic cylinders 
and the resonator enclosing it, yielding considerably di- 
minished widths of the resonances. Still, in the vicinity of 
the Dirac frequency the density of the resonances is low 
such that they can be resolved. In fact, there the local 
density of states tends to zero linearly with the distance 
I / — /d I of the excitation frequency / from the Dirac fre- 
quency fu^d&dL. Note that only in this frequency range 
the properties of the microwave Dirac billiard are well 
described by the Dirac equation, that is, only there it 
can simulate a Dirac billiard. The associated frequency 
range has been determined in Ref. 17 using the relation 
between the local density of states and reflection spec- 
tra measured with one antenna placed in the interior of 
the photonic crystal. It extends from / ~ 13.3 GHz 
to / ~ 14.3 GHz. In this range of separated resonances 
the corresponding wave functions could be measured and 
their properties, such as the occurence of edge states pre- 
dicted for graphene billiards 24 , were investigated. 



The wave functions of a quantum billiard are directly 
related to the electric field strength distributions in the 
corresponding microwave billiard at the resonance fre- 
quencies^. Thus, their intensities can be determined ex- 
perimentally by measuring the electric field intensity dis- 
tributions. This was done with the help of the so-called 
perturbation body metho d 32 i 43 i 44 , which is based on the 
Maier- Slater theorem 4 ^. A small perturbation body is 
inserted into the cavity. This induces a shift of the reso- 
nance frequencies the size of which depends on the elec- 
tric and magnetic field strengths at the position of the 
perturbation body. Since we are only interested in the 
intensity distribution of the wave functions of the cor- 
responding Dirac billiard, we used a perturbation body 
made from magnetic rubber—. In order to obtain the 
spatial distribution of the electric field strength the per- 
turbation body was moved in the closed cavity in 2 mm 
steps between the 13 rows of cylinders which are paral- 
lel to the long sides of the billiard. This was done with 
an external guiding magnet, which was placed on a com- 
puter controlled positioning unit. For each position the 
transmission between the two antennas marked in Fig. [71 
was measured. 

Figure [9] shows exemplarily ten electric field intensity 
distributions measured at resonance frequencies around 
the Dirac frequency. The labels of the panels correspond 
to the associated resonances in Fig. [8j The intensity pat- 
terns shown in panels a and j are extended over the whole 
area of the microwave billiard. The structure of these in- 
tensity distributions resembles that of modes of the cor- 
responding empty rectangular quantum billiard with one 
excitation parallel to the short billiard edge and several 
parallel to the long one. They are similar to the wave 
functions shown in Fig. 3 (b) of Ref. 27 . The intensity 
distributions shown in panels c - i are localized close to 
a long edge of the rectangular billiard, i.e., there the 
intensity is highest and it is much smaller in the inte- 
rior, especially in panels e - h. At the long sides the 
photonic crystal terminates with straight edges, at the 
short sides with zigzag edges. These are close to, re- 
spectively, the zigzag edges and the armchair edges in 
the corresponding graphene billiard defined by the void 
structure of the photonic crystal. An example is shown 
in Fig. [TJ where the two sublattices are depicted as, re- 
spectively, red and blue dots. States which are localized 
near a short edge of the microwave billiard were not ob- 
served 16 , whereas altogether 16 corresponding to a se- 
quence of resonances around the Dirac frequency were 
found to be localized near a long edge. This is in good 
agreement with the results for the associated rectangular 
graphene billiard 2 ^. There, 14 consecutive states around 
the Dirac point, which are all localized along a zigzag 
edge, so-called edge states, are predicted. Note that the 
boundary conditions of the microwave billiard and the 
corresponding graphene billiard differ. For the latter the 
wave functions vanish on one sublattice at a zigzag edge, 
and on both at an armchair edge, whereas for the former 
they are zero along the side walls, which are arranged 
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FIG. 9: (Color online) The electric field intensity distribu- 
tions measured at the resonance frequencies in Fig. [8] These 
are marked with arrows and the labels ranging from a— Kj. 
The blue (darkest) color corresponds to the lowest, the red 
(brightest) to the highest intensity. 

slightly apart from the edges of the graphene billiard (see 
Fig. [I]). In spite of these differences both systems exhibit 
edge states corresponding to consecutive resonance fre- 
quencies around the Dirac frequency, which all are local- 
ized near a zigzag edge. They remind of the edge states 



observed experiment ally 28 r— in graphene flakes, where 
electrons are confined by a potential barrier which, in 
distinction to that along the boundary of the Dirac bil- 
liard, is of finite height, and recently also in a magnetic 
photonic crystal 48 . 

III. SUMMARY 

In summary, we have observed experimentally the 
anomalous 1/L scaling behavior of the transmission of 
electromagnetic waves through a photonic crystal, which 
were excited at the Dirac frequency and sent onto the 
crystal in directions close to the TK direction, that is, 
near a conically shaped singularity of the photonic band 
structure. Furthermore, for the first time we imple- 
mented an experimental setup representing a microwave 
analog of a Dirac billiard, that is a relativistic quantum 
billiard. We have observed experimentally in the vicin- 
ity of the Dirac point states localized at edges of the 
Dirac billiard, which are similar to the edge states occur- 
ring at zigzag edges in graphene flakes. These similarities 
between the properties of waves propagating through a 
photonic crystal and those of graphene corroborate ex- 
perimentally that not only the linear dispersion relation 
but also the unusual properties of graphene associated 
with it are entirely due to the symmetry properties of its 
lattice. Furthermore, this demonstrates that microwave 
billiards containing a photonic crystal provide an exper- 
imental setup for the investigation of properties of Dirac 
billiards of different shapes and with varying boundary 
conditions. These systems are particularly suited for the 
investigation of the spectral properties of Dirac billiards 
since the measurements can be performed with super- 
conducting microwave billiards at liquid helium temper- 
ature. The aim is the study of the effect of the shape 
of a Dirac billiard and the edge structure of the enclosed 
photonic crystal on its spectral properties. Another inter- 
esting aspect concerning properties of the wave functions 
is the scarring phenomenon which has been seen numer- 
ically in graphene flakes with the shape of a stadium 47 . 
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